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摘要 HODHA (Rhododendron delavayi) 花 开 花 至 凋谢 进程 中 的 代谢 产物 差异 及 其 


通路 ， 该 文采 用 LC-MS/MS 技术 对 


花苞 期 、 开 裂 期 、 传 粉 期 、 盛 开 期 、 衰 老 期 和 凋谢 期 


的 化 学 成 分 进 


行 非 靶 向 代谢 组 学 分 析 。 结 果 表明 : d) REER 973 种 代谢 物 ， 主 要 包含 


黄酮 类 、 


E 


Fc < 0.5) 


dg or 
者 员 集 


物 合成 。 


WAAR (L-phenylalanine) 、 反 式 肉 桂 酸 (trans-cinnamic acid) ~ Æ 
素 Cnaringenin) 、 对 香 豆 酰 基 莽 草酸 Cp-coumaroyl shikimic acid) ~ 


明 样 本 间 代 谢 物 存在 差异 
变量 分 忆 


* EJ, 其 中 花色 期 至 开裂 期 差异 代谢 物 的 表达 主要 呈现 下 调 , 而 进入 衰老 期 和 凋谢 期 后 
差异 代谢 物 的 表达 主要 呈 
(P < 0.01) 通路 3 条， 包括 茶 丙 素 类 生物 合成 、 植 物 激素 的 生物 合成 和 类 黄酮 生 


机 酸 、 酚 酸 类 


、 氨 基 酸 及 其 衍生 物 、 脂 类 、 生 物 碱 等 。(2) 主 成 分 分 析 (PCA) 
， 结 合 正 交 偏 最 小 二 乘 判 别 分 析 (OPLS-DA) 、t 检验 的 P 值 和 


的 差异 倍数 (fold-change, Fe) 筛选 差异 代谢 物 (VIP > 1，P<0.05，Fc >2 或 
， 涉 及 591 种 ， 在 马 绢 杜 鹏 花期 进入 衰老 期 和 凋谢 期 后 差异 代谢 物 数量 和 表达 


FH 


WEW. (3) KEGG 注释 到 68 条 代谢 通路 ， 其 中 差异 代谢 物 极 显 


分 生 


(4) 结合 葵 丙 素 类 、 黄 酮 类 等 有 效 成 为 合成 通路 共 篇 选 到 10 种 代谢 物 包括 茶 
耳 酮 Cchalcone) 、 柚 皮 


阿 魏 酸 (ferulic acid) ~ 


松柏 醇 (coniferyl alcohol) 、 芥 子 酸 Csinapic acid) . 4€ J & ff Csyringin) . WHA (quercetin). 


IEAk, AARDT E] 28 RNC] SAT RREA RR E s RARE BS TE DIR Z2 EE AAA, 
而 黄酮 类 化 合 物 生物 合成 逐渐 减弱 , Xt EE ORE TT 26 jr ECCE RT BEOGE s RERBA EE 
的 调控 作用 。 该 研究 为 马 继 杜 鹏 花 开 花 至 凋谢 进程 中 的 有 效 成 分 代谢 途径 活性 物质 的 研究 提 
供 了 代谢 组 学 基础 ， 为 进一步 研究 马 织 杜 鹏 花花 期 调控 的 分 子 机 理 提供 参考 。 
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Abstract: In order to analyze the differences of metabolites and their pathways from flower 
opening to withering of Rhododendron delavayi, we used LC-MS/MS material separation and 
identification technique to non-targeted the chemical components of bud stage, dehiscence stage, 
pollination stage, blooming stage, senescence stage and withering stage. The results were as 
follows: (1) A total of 973 kinds of metabolites were detected, mainly including flavonoids, 
organic acids, phenolic acids, amino acids and their derivatives, lipids, alkaloids and so on. (2) 
Principal component analysis (PCA) shows that there were differences in metabolites among 
samples. Combined with orthogonal partial least squares discriminant analysis (OPLS-DA), P 
value of t-test and fold change (Fc) of univariate analysis, differential metabolites (VIP > 1, P < 
0.05, Fc > 2 or Fc < 0.5) were screened out, involving 591 species, the quantity and expression of 
differential metabolites increased significantly after the Rhododendron delavayi flower stage 
entered the senescence stage and the withering stage, in which the expression of differential 
metabolites from bud stage to dehiscence stage was mainly down-regulated, while the expression 
of differential metabolites was mainly up-regulated after entering senescence stage and withering 
stage. (3) 68 differential metabolic pathways were annotated by KEGG, of which 3 pathways were 
significantly enriched with differential metabolites (P < 0.01), including phenylpropanoids 
biosynthesis, plant hormone biosynthesis and flavonoid biosynthesis. (4) Based on the analysis of 
biosynthetic pathways of phenylpropanoids, flavonoids and other effective components, 10 kinds 
of differential metabolites were screened from flowering to withering of Rhododendron delavayi, 
including L-phenylalanine, trans-cinnamic acid, chalcone, naringenin, p-coumaroyl shikimic acid, 
ferulic acid, coniferyl alcohol, sinapic acid, syringin and quercetin. In addition, the differential 
metabolites of effective components showed that phenylpropanoids biosynthesis and metabolism 
activities gradually increased with the development of Rhododendron delavayi, while flavonoid 
biosynthesis gradually decreased. These key differential metabolites may play an important role in 
regulating the development of Rhododendron delavayi. This research provides a metabonomic 
basis for the study of effective components in the metabolic pathway during the flowering and 
withering process of Rhododendron delavayi, and provides a reference for the further study of the 
molecular mechanism of flowering regulation of Rhododendron delavayi. 

Key words: Rhododendron delavayi, flowering period, metabolome, LC-MS/MS, metabolic 
pathway 


IAES (Rhododendron delavayi) 属于 常 绿 灌木 或 小 乔木 ， 因 其 花色 鲜艳 、 形 状 酷似 
ILKKA, 主要 分 布 于 西南 地 区 , 包括 四 川西 南部 、 云 南 全 省 和 贵州 西北 部 等 ( 张 
K 8$, 20150 ， 其 中 贵州 西北 部 以 百 里 杜 鹏 的 普 底 和 人 金 坡 景区 分 布 最 为 广泛 ， 为 主要 的 
建 群 种 和 优势 种 〈 李 甘 洁 和 陈 训 ，2005) 。 每 年 3 月 中 旬 前 完成 初 花期 ( 花 泊 开始 开张 到 
10% 人 花序 开放 ) ; 3 月 下 旬 到 4 月 初 进 入 为 盛开 期 (10% 花序 开放 到 10% 花 序 出 现 凋谢 ) ; 4 
J E d P 709946 FF tH BLUR EO , S ZEREBS SEA EHI 25 一 
30 天 。 研 究 马 继 杜 鹏 花 开 花 至 凋谢 进程 中 的 代谢 组 分 变化 ， 对 延长 和 调节 马 继 杜 鹏 花期 ， 
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$54] AT 3 PEEEBS AN AEA E IAD AN AR SRRA BURGESEBS 
研究 意义 。 目 前 关于 马 继 杜 鹏 花 主 要 包括 花 特 征 和 花粉 性 状 的 研究 ， 以 百 里 杜 岗 保护 区 马 继 
社 鹏 为 实验 材料 研究 花 特 征 分 化 和 资 众 源 分 布 情况 〈 胡 德 美 等 ，2021) ; 用 显微镜 观察 马 继 
杜 鹏 花粉 形态 大 小 、 外 壁 纹饰 、 孔 沟 形 态 等 性 状 〈 毛 加 宁 ，2000) 。 同 时 ， 在 马 继 杜 鹏 花花 
期 调控 的 研究 则 主要 集中 于 营养 元 素 、 植 物 激素 和 温度 处 理 等 方法 ， 包 括 用 不 同 浓度 的 N、 
P, K, Ca 等 营养 元 素 处 理 马 绢 杜鹃 观察 对 花芽 发 育 和 花期 的 影响 ( 洪 鳗 等 ，2010) ; 在 
马 绢 杜 鹏 花芽 发 育 初 期 和 开花 前 期 喷 施 外 源 激素 对 始 花 期 、 盛开 期 、 末 花期 和 开花 数 的 影响 
(RRR 等 ，2010) ;通过 不 同 浓度 的 生长 调节 物质 及 矿物 质 营养 处 理 马 继 杜 鹏 的 花蕾 ， 
探索 马 继 杜 鹏 花花 期 调控 方法 〈 策 保 超 ，2015) 。 目 前 仅 有 少量 关于 马 织 杜 鹏 茎 和 根 代谢 物 
的 研究 报道 ， 利 用 各 种 色谱 技术 从 马 继 杜 鹏 的 茎 分 离 鉴 定 得 到 10 种 化 合 物 〈 徐 金 金 等 ， 
2012) ， 然 而 马 织 杜 鹏 花 的 发 育 进 程 中 代谢 物 和 代谢 物 途径 研究 还 未 见 相 关 的 报道 。 
近年 来 , 植物 代谢 组 学 研究 已 经 被 广泛 地 应 用 于 分 析 果 树 生 长 、 发 育 、 果 实 成 熟 过 程 中 
的 代谢 物 变 化 (Aaron et al., 2008; Zhang et al., 2011; Li etal., 2020) ,目前 也 开始 应 用 于 而 
究 植 物 花 形成 、 发 育 的 代谢 机 制 研究 ( 贾 岩 等 ，2017; DLA 等 ，2020) ， 采 用 气相 色谱 
-质谱 法 技术 分 析 款 冬 花 不 同 发 育 阶段 代谢 组 ， 发 现款 冬 花 董 发 育 过 程 中 部 分 特征 代谢 物 呈 
~ 现 显著 的 变化 〈 薛 水 玉 等 ，2012〉; 利用 转录 组 和 代谢 组 研究 铁皮 石 解 花 的 发 育 机 制 ， 发 
现 次 生 代谢 和 碳水 化 合 物 的 变化 与 花 的 发 育 有 关 (He et al., 2020) ;采用 GC-MS 技术 检测 
枇杷 花 发 育 进 程 中 的 代谢 物质 , 揭示 影响 花 发 育 进程 的 生理 与 代谢 机 制 ( 徐 红 霞 等 , 2020)。 
代谢 组 学 对 植物 花 的 发 育 代谢 产物 及 代谢 通路 研究 提供 了 方法 。 

本 研究 以 贵州 省 毕节 市 大 方 县 百 里 杜 鹏 风景 区 为 研究 区 域 ， 以 马 织 杜 鹃 花 为 研究 对 象 ， 
选取 花苞 期 、 开 裂 期 、 传粉 期 、 盛开 期 、 衰 老 期 和 凋谢 期 为 材料 ， 采 用 非 革 向 代谢 组 学 
(LC-MS/MS) 分 析 方 法 ,通过 马 绢 杜 鹏 花 开花 至 凋谢 进程 中 6 个 阶段 的 代谢 产物 及 代谢 通 
路 研究 ， 拟 探讨 以 下 问题 : CD 马 继 杜 鹏 花 开 花 至 凋谢 进程 中 不 同 阶段 代谢 物 并 且 是 如 何 
变化 的 ? (20 分 析 马 绢 杜 鹏 花 差 异 显 著 的 代谢 物 和 相关 的 代谢 通路 ; (3) Ox h E S HR 
赐 花 发 育 相关 的 重要 代谢 物 ， 该 研究 结果 为 马 继 杜 鹏 花花 期 调控 提供 理论 基础 和 科学 依据 。 


1 材料 与 方法 


E L1 试验 材料 
"e 试验 材料 为 贵州 省 毕节 市 大 方 县 百 里 杜 鹏 风景 区 马 继 杜 鹏 (Rhododendron delavayi 
€) Franch.) . 


2020 年 选取 长 势 基本 一 — EA A PUE 
传粉 期 (3 月 30 日 ) 、 盛 开 期 (4 月 衰老 期 (4 月 10 日 ) 、 调 谢 期 (4 月 15 日 ) ， 
从 6 株 树 上 采集 适量 马 继 杜 鹏 花样 本 ， o 回 实验 室 进行 代谢 组 学 检测 , 每 个 时 
期 取 6 个 生物 学 重复 样本 。 

1.2 马 继 杜鹃 花 代谢 物 提 取 

基于 液 质 联 用 (LC-MS/MS) 技术 (Warwick et al., 2011; Want et al., 2010) 3t£T 4ESU [n] 
代谢 组 学 研究 。 取 100 mg 液 氮 研 磨 的 马 织 杜 鹏 花样 本 ， 置 于 EP 管 中 ， 加 入 500 uL 0.1% 
甲酸 的 80% 甲 醇 水 溶液 ， 涡 旋 震 荡 ， 冰 浴 静 置 5 min，15 000 rpm、4 CHÙ 10 min， 取 100 
uL 的 上 清 加 超 纯 去 离子 水 稀释 至 甲醇 含量 为 $3%, 并 置 于 离心 管 中 15 000 g、4 "C Aù 10 min, 
收集 上 清 液 ， 进 样 LC-MS/MS 进行 检测 分 析 。 

从 每 个 实验 样本 中 取 等 体积 样本 混 匀 作为 QC 样本 。 空 白 样本 为 含 0.1% 甲 酸 的 5396 FH 
醇 水 溶液 代替 实验 样本 ， 前 处 理 过 程 与 实验 样本 相同 。 

1.3 非 靶 向 代谢 组 学 (LC-MS/MS) 检测 

采用 Thermo Hyperil Gold column 色谱 柱 〈Cig) ， 柱 温 40 C, FUR 0.2mL/min. ES 

子 模式 : 流动 相 A: 0.1% 甲 酸 ， 流 动 相 B: 甲醇 ; 负离子 模式 : 流动 相 A: SmM MRE (pH 
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9.00 ， 流 动 相 B: 甲醇 ， 梯 度 洗 脱 (0 — 1.5 min，98%A，2%B; 1.5 ~ 12min, 10094B; 
12 ~ 14min, 100%B; 14 ~ 14.1 min, 9896A, 296B; 14.1 ~ 17.0min，98%A，2%B ) 。 
进 样 量 为 10 uL. AH Thermo Vanquish UHPLC 液 相 色 谱 Thermo QE 系列 质谱 分 别 进行 正 、 
负离子 模式 采集 。 扫描 范围 选择 m/z 70 ~ 1050; ESI 源 的 设置 如 下 : Spray Voltage: 3.2 kV; 
Sheath gas flow rate: 35arb; Aux Gasflow rate: 10arb; Capillary Temp: 320 °C. Polarity: 
positive; negative; MS/MS 二 级 扫描 为 data-dependent scans. f$ 10 个 分 析 样 本 进行 一 次 QC 
样本 ， 评 价 整个 实验 过 程 中 仪器 状态 的 稳定 性 和 检测 数据 结果 的 准确 性 。 

1.4 数据 分 析 

将 采集 到 的 原始 数据 导入 CompoundDiscoverer3.1 (CD) 软件 中 ， 通 过 保留 时 间 、 质 荷 
比 对 数据 进行 过 滤 ， 然 后 根据 保留 时 间 偏 差 0.2 min 和 质量 偏差 5 ppm 进行 峰 对 齐 ， 随 后 根 
据 质 量 偏 差 5 ppm、 信 号 强度 偏差 30% (Dai etal., 2017) 等 信息 进行 峰 提 取 ， 同 时 对 峰 面积 
进行 定量 ， 再 整合 目标 离子 ， 然 后 通过 分 子 离子 峰 和 碎片 离子 进行 分 子 式 的 预测 并 与 
mzCloud. mzVault 和 MassList 数据 库 进 行 比 对 ， 用 blank 样本 去 除 背 景 离子 ， 并 对 定量 结 
果 进 行 归 一 化 ， 最 后 得 到 数据 的 鉴定 和 定量 结果 。 

使 用 MetaX (Wen et al., 2017) 软件 对 数据 进行 对 数 转换 处 理 ， 将 所 有 实验 样本 和 QC 
样本 导入 SIMCA14.1 软件 进行 多 元 统计 分 析 《 阿 基业 ，2010) ， 采 用 主 成 分 分 析 (PCA) 
观察 各 样本 的 整个 过 程 的 稳定 性 ， 然 后 进一步 使 用 有 监督 的 正 交 偏 最 小 二 乘 判别 分 析 
COPLS-DA) 来 区 分 各 组 间 代 谢 物 的 差异 性 。 根 据 多 变量 OPLS-DA 模型 分 析 中 变量 权重 
(VIP>1) 、 结 合 上 检验 的 已 值 CP < 0.050 和 单 变量 分 析 差 异 倍数 CFold-change, Fc >2 
或 Fc <0.5) 筛 选 差 异 代谢 物 。 维 恩 图 由 TB tools 0.6 软件 完成 ,代谢 通路 富 集 分 析 基 于 KEGG 
数据 库 。 


2 结果 与 分 析 


2.1 马 继 杜 鹃 花 开花 至 凋谢 过 程 的 形态 变化 

将 马 继 杜 鹃 花 开 花 至 凋谢 过 程 分 为 6 个 时 期 ， 如 图 1 所 示 。 花 区 期 (A) : ERF, K 
现 为 花苞 状态 ， 内 部 花药 还 未 成 熟 ， 雄 营 未 开裂 开裂 期 B): 花苞 已 经 初步 打开 ,雄蕊 
刚 开 有 裂 ， 雌 营 柱 头 上 未 有 花粉 ， 传粉 期 OOO : 花苞 较为 开发 ， 雄 蕉 完 全 裂 开 ， 雌 营 柱 头 上 
已 经 结合 覆盖 花粉 ， 盛开 期 CODO : 花 东 呈现 完全 盛开 ， 颜 色 艳 丽 ， 雄 蕉 已 经 完全 开裂 ， 且 
花粉 已 经 全 部 裂 出 ， 雌 营 柱 头 上 花粉 已 经 不 多 ， 传 粉 受精 已 完成 ;衰老 期 CE) : 花 末 开始 
凋谢 ， 出 现 蓉 芒 状态 ， 调 谢 期 (F) : REES, AME. 
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A. 花苞 期 ，B. 开裂 期 ，C. 传粉 期 ，D. 盛开 期 ，E. 衰老 期 ，F. 凋谢 期 。 


A. Bud stage; B. Dehiscence stage; C. Pollination stage; D. Blooming stage; E. Senescence stage; F. Withering 


stage. 
图 1 马 继 杜 鹏 花 开花 至 凋谢 过 程 的 形态 变化 

Fig.l Morphological changes of Rhododendron delavayi from flowering to withering 
2.2 马 继 杜 鹃 花样 本 主 成 分 分 析 

为 分 析 马 继 杜 鹏 花 各 个 时 期 间 代 谢 组 的 相互 关系 ， 采 用 主 成 分 分 析 CPCAO 通过 建立 数 
学 模型 进行 正 交 变换 简化 降 维 分 析 样 本 间 代 谢 谱 的 内 部 结构 关系 。 马 继 杜 鹏 花 主 成 分 分 析 
(图 2: A 表示 正 离子 模式 ，B 表示 负离子 模式 ) ， 结 果 表明 不 同时 期 代谢 物 各 组 间 存 在 明 
显 的 分 离 趋势 ， 说 明 每 个 样本 的 数据 处 理 结果 是 可 信 的 ， 每 个 样本 之 间 存 在 明显 差异 ， 且 组 


了 


内 生物 学 重复 重合 效果 好 , 表明 在 实验 过 程 中 数据 组 的 重 现 性 好 , 可 以 用 于 后 续 差 异 代谢 物 
分 析 。 花 区 期 和 开裂 期 、 传 粉 期 、 盛 开 期 代谢 物 的 变化 相对 较 大 ， 衰 老 期 与 凋谢 期 的 代谢 物 
相对 最 大 。 


QC. 质 控 样本 ，B. 花苞 期 ，PD. 开裂 期 ，P. 传粉 期 ，AP. 盛开 期 ，S. 衰老 期 ; F. 凋谢 期 。 下 同 。 
QC. Quality control sample; B. Bud stage; PD. Dehiscence stage; P. Pollination stage; AP. Blooming stage; S. 


Senescence stage; F. Withering stage. The same below. 

图 2 马 织 杜 鹏 花 开花 至 凋谢 进程 中 代谢 物 PCA 模型 得 分 图 

Fig.2 PCA model scores of metabolites from flowering to withering of Rhododendron delavayi 
为 筛选 出 马 继 杜 鹏 开花 至 凋谢 进程 中 的 差异 代谢 物 ， 采 用 正 交 偏 最 小 二 乘 判别 分 析 

COPLS-DA) 区 分 6 个 时 期 代谢 物 的 总 体 差异 ， 结 果 如 图 3 所 示 。 马 继 杜 鹏 开花 至 凋谢 进 

程 中 ， 花 奋 期 、 开 和 裂 期 、 传 粉 期 、 盛 开 期 、 衰 老 期 、 凋 谢 期 6 个 时 期 的 代谢 物 明显 分 离 ， 表 

明和 名 组 样本 代谢 物 存 在 显著 性 的 差异 变化 。 该 模型 质量 参数 R? 和 QQ 均 大 于 0.9， 说 明 当 前 

OPLS-DA 分 析 模 型 稳定 可 靠 ， 具 有 良好 的 预测 能 力 ， 有 利于 寻找 潜在 的 马 绢 杜 鹏 开花 至 凋 

谢 进程 中 的 差异 代谢 物 。 
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A. 正 离子 模式 ，B. 负离子 模式 。 
A. Positive ion mode; B. Negative ion mode. 
图 3 BARBERER VATER P NI OPLS-DA 得 分 图 
Fig.3 OPLS-DA scores of metabolite from flowering to withering of Rhododendron 


delavayi 

22 BAHR ER HHN AA 

E R 3 ARSTER 46 8 UAE E ER PARARE, RAH LC-MS/MS M 5 2E ERSTE 6 
个 时 期 的 代谢 物 多 峰 图 ， 花 奋 期 代谢 物色 谱 图 见 图 4， 其 它 时 期 代谢 物 峰 图 见 附件 1。 共 检 
测 到 973 种 代谢 物 ， 正 离子 模式 587 种 ， 负 离子 模式 386 种 ， 其 中 黄酮 类 114 种 ， 有 机 酸 
114 种 ， 酚 酸 类 112 种 ;氨基 酸 及 其 衍生 物 75 种 ， 脂 类 61 b. EWT STU W 45 种 ; 
核 背 酸 及 其 衍生 物 43 种 ; 糖 类 及 其 衍生 物 36 种 ， 茶 丙 素 类 33 fb; 醇和 多 元 醇 33 RB ERA 
11 种 ;省 体 3 种， 其 他 类 236 种 。 
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A. 正 离子 模式 ，B. 负离子 模式 。 
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A. Positive ion mode; B. Negative ion mode. 
图 4 LC-MS/MS 代谢 物 检 测 多 峰 图 
Fig.4 LC-MS/MS metabolite detection multimodal graph 

2.3 差异 代谢 物 分 析 

根据 OPLS-DA 模型 的 变量 重要 性 投影 值 第 选 马 绢 杜鹃 花 6 个 时 期 的 差异 代谢 物 (VIP > 
1、P<0.05 和 Fc (fold change) >2 或 fold change « 0.5) 。 从 花苞 期 、 开裂 期 、 传 粉 期 、 
盛开 期 、 衰 老 期 到 凋谢 期 进行 相 邻 两 个 时 期 的 代谢 组 比较 ,结果 发 现 : 马 继 杜 鹃 花 开花 到 凋 
谢 进程 中 共有 591 种 差异 代谢 物 。 其 中 , 花苞 期 与 开裂 期 有 170 种 差异 代谢 物 ( 正 离子 模式 
87、 负 离子 模式 83) ，48 种 上 调 、122 PE P Jd; 开裂 期 与 传粉 期 有 144 种 差异 代谢 物 〔 正 
离子 模式 89、 负 离子 模式 55) ，68 种 上 调 、76 种 下 调 ; 传粉 期 与 盛开 期 有 111 种 差异 代谢 
物 〈 正 离子 模式 67、 负 离子 模式 44) ，75 种 上 调 、36 种 下 调 ;盛开 期 与 衰老 期 有 284 种 差 
异 代 谢 物 〈 正 离子 模式 184、 负 离子 模式 100) ，177 种 上 调 、107 种 下 调 ; 衰老 期 与 凋谢 期 
1314 种 差异 代谢 物 〈 正 离子 模式 180、 负 离子 模式 134) ，260 种 上 调 、54 种 下 调 。 在 马 
继 杜 鹃 花期 进入 衰老 期 和 凋谢 期 后 差异 代谢 物 数量 明显 增多 于 开花 过 程 , 且 大 部 分 是 表达 量 
显 普 上 升 的 差异 代谢 物 。 而 在 马 继 杜 赐 花苞 期 至 开裂 期 过 程 中 , 差异 代谢 物 中 大 部 分 的 表达 
量 显 著 下 调 。 差 异 代谢 物 维 恩 图 分 析 表 明 (图 5) : 绝 大 部 分 (246.8 %) 的 差异 代谢 物 都 
只 出 现在 一 个 相 邻 变化 时 期 没有 任何 差异 代谢 物 出 现在 整个 开发 至 凋谢 过 程 。 


A. 上 调 代 谢 物 ，B. 下 调 代 谢 物 。 
A. Up-regulated metabolite; B. Down-regulated metabolite. 
图 5 马 继 杜 鹏 花 开花 至 凋谢 进程 中 差异 代谢 物 维 恩 图 


Fig. 5 Venn diagram of differential metabolites from flowering to withering of 


Rhododendron delavayi 

2.4 代谢 通路 分 析 

通过 KEGG Pathway 富 集 分 析 显 著 差 异 代 谢 物 CKanehisa & Goto, 2000; Jun Rao et al., 
2014; Lin Hong et al., 2014) ， 如 图 6 所 示 。 注释 到 68 条 代谢 通路 , 显著 代谢 通路 (P< 0.05) 
要 有 6 条 ， 分 别 是 茶 丙 素 类 生物 合成 ， 植 物 激素 的 生物 合成 ， 类 黄酮 生物 合成 ， 精 氨 酸 和 及 
氨 酸 代谢 ， 茶 丙 氨 酸 、 栈 氨 酸 和 人 色 和 氨 酸 生物 合成 鸟 氨 酸 、 赖 氨 酸 和 烟 酸 生 物 碱 的 生物 合成 。 
其 中 ， 极 显著 的 代谢 通路 〈P< 0.01) 包括 茶 丙 素 类 生物 合成 、 植 物 激素 的 生物 合成 和 类 黄 
酮 生物 合成 。 
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HERALA VERUS ERER EA k Biosynthesis of terpenoids and steroids : 
不 饱和 脂肪 酸 的 生物 合成 Biosynthesis of unsaturated fatty acids ; 
氨 素 代谢 Nitrogen metabolism | 
生物 素 代谢 Biotin metabolism | 
泛酸 盐 和 CoA 生 物 合成 Pantothenate and CoA biosynthesis | 
FNRA EHDA EMA R Biosynthesis of alkaloids derived from shikimate pathway | 。 
叶酸 碳 库 One carbon pool by folate ，。 
淀粉 和 蔗糖 代谢 Starch and sucrose metabolism ; 
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氨 酰 tRNA 生物 合成 Aminoacyl-tRNA biosynthesis + 


EAR} Purine metabolism : 

甘氨酸 、 丝 氨 酸 和 苏 氨 酸 代谢 Glycine, serine and threonine metabolism * 

黄酮 和 黄酮 醇 生物 合成 Flavone and flavonol biosynthesis ^ 

鸟 氨 酸 、 赖 氨 酸 和 烟 酸 生物 碱 的 生物 合成 Biosynthesis of alkaloids derived from ornithine, lysine and nicotinic acid | 
葵 两 氨 酸 、 栈 氨 酸 和 色 氨 酸 生 物 合成 Phenylalanine, tyrosine and tryptophan biosynthesis | 

精 氨 酸 和 及 氨 酸 代谢 Arginine and proline metabolism ] 

类 黄酮 生物 合成 Flavonoid biosynthesis | 

植物 激素 的 生物 合成 Biosynthesis of plant hormones | 

茶 丙 素 的 生物 合成 Biosynthesis of phenylpropanoids . 


Fig.6 Enrichment analysis of main metabolic pathways 
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Klo 马 继 杜 鹏 花 开 花 至 凋谢 进程 中 差异 代谢 物 主要 代谢 通路 富 集 分 析 


of differential metabolites from 


flowering to withering of Rhododendron delavayi 


本 研究 根据 代谢 组 学 分 析 差 异 代谢 物 与 代谢 通路 的 关系 


, 挖掘 与 马 继 杜 赐 花 开花 至 凋谢 


相关 的 差异 代谢 物 ， 推 测 代谢 物 可 能 的 作用 机 制 。 结 合 代谢 通路 图 〈 图 7) ， 共 筛选 到 10 


HERRA, GARNAR (L-phenylalanine) 、 反 式 肉桂 酸 (trans-cinnamic acid) 、 查 


耳 酮 (chalcone) 、 柚 皮 素 (naringenin) 、 对 香 豆 酰 基 莽 草酸 (p-coumaroyl shikimic acid) 、 
阿 魏 酸 (ferulic acid)、 松 柏 醇 (coniferyl alcohol)、 芥 子 酸 (sinapic acid)、 紫 丁香 苷 (Csyringin )、 
HZR Cquercetin) . HAIRS 6 个 时 期 ， 代 谢 物 含量 表达 如 图 6 所 示 ， 茶 丙 素 类 生 
物 合成 相关 的 差异 代谢 物 包 括 茶 丙 氨 酸 (L-phenylalanine)、 反 式 肉桂 酸 (trans-cinnamic acid)、 
对 香 豆 酰 莽 草酸 (p-coumaroyl shikimic acid) ~ TA4HB£ Cconiferyl alcohol) ~ FFIR Csinapic 


acid) 、 紫 丁香 苷 〈syringin) 在 衰老 期 (8S) 和 凋谢 期 CF) 


两 个 时 期 均 为 上 调 表达 ， 特 别 


是 在 凋谢 期 P 达到 最 高 值 ， 而 黄酮 类 生物 合成 相关 的 代谢 物 柚 皮 素 Cnaringenin) . H 


素 (quercetin) 和 查 尔 酮 (chalcone) 随 着 花 的 生长 发 育 含量 


季 果 冬 芽 分 化 关键 时 期 ，CCC 处 理会 诱导 茶 丙 素 合成 和 类 黄 


呈现 下 调 表 达 。 在 ' 夏 黑 ' 葡 萄 二 
酮 合成 通路 中 部 分 基因 的 上 调 


或 下 调 表 达 , 通过 改变 与 成 花 相 关 的 基因 表达 和 次 级 代谢 产物 的 合成 和 代谢 来 实现 促进 成 花 


(时 晓 芳 等 ，2021) 。 转 录 组 学 分 析 不 同 发 育 阶段 款 冬 叶 ， 
丙 素 类 生物 合成 相关 基因 表达 量 不 断 增 多 《县 佳 慧 等 ，201 


发 现 叶 片 生长 发 育 的 过 程 中 茶 
8) 。 马 继 杜 鹏 花 开 花 至 代谢 进 


程 中 葵 丙 素 类 代谢 物 表 达 量 逐渐 升 高， 而 黄酮 类 代谢 物 表 达 量 逐渐 降低 , 葵 丙 素 类 和 黄酮 类 


代谢 物 表 达 可 能 与 马 继 杜 鹏 花 发 育 有 关 。 
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图 7 茶 丙 素 类 、 黄 酮 类 等 有 效 成 分 生物 合成 通路 


Fig.7 Biosynthesis pathway of active components about phenylpropanes and flavones 
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3 讨论 与 结论 


3.1 马 继 杜 鹏 花 差 异 代谢 物 

代谢 物 是 生物 体 表 型 的 基础 ， 能 帮助 我 们 更 直观 有 效 地 了 解 生物 学 过 程 及 其 机 理 〈 熊 欢 
欢 等 ，2019) 。 采 用 GC-MS 技术 对 枇杷 花 整 个 发 育 进程 中 代谢 物质 分 析 ，C/N 含量 比 是 决 
定 花 发 育 进程 的 重要 因素 之 一 〈 徐 红 霞 等 ，2020) 。 长 期 以 来 马 继 杜 鹏 花花 期 调控 一 直 都 
是 研究 热点 , 马 继 杜 鹏 花 的 开花 至 凋谢 进程 是 十 分 复杂 的 生理 过 程 , 该 过 程 可 能 与 代谢 物种 
类 和 含量 密切 相关 。 本 研究 采用 非 靶 向 代谢 组 学 方法 , 研究 马 绢 杜 鹏 花 开花 至 凋谢 的 6 个 时 
期 〈 花 区 期; 开裂 期 , 传粉 期 盛开 期 ;衰老 期 和 凋谢 期 ) 代谢 物 变 化 。 通 过 代谢 组 学 分 析 ， 
共 筛 选 获得 591 种 差异 表达 代谢 物 , 主要 包括 黄酮 类 、 有 机 酸 、 酚 酸 类 、 氨基 酸 及 其 衍生 物 、 
脂 质 和 生物 碱 等 , 这 些 差异 代谢 物 在 花 的 生长 发 育 不 同时 期 表达 情况 不 相同 ,， 随 着 马 继 杜 鹏 
花 的 生长 发 育 差异 代谢 物 在 花苞 期 至 盛开 期 逐渐 减少 , 进入 衰老 期 和 凋谢 期 后 显著 增加 ， 而 
上 调 代谢 物 在 整个 开花 至 凋谢 进程 中 逐渐 增加 , 并 且 在 凋谢 期 达到 最 大 值 , 下调 代谢 物 在 花 
区 期 至 怀 开 期 逐渐 减少 ， 衰老 期 增加 ， 凋谢 期 减少 ， 本 研究 马 继 杜 鹃 花 的 不 同 发 育 阶段 差异 
代谢 物种 类 和 表达 可 能 与 花 的 开花 至 凋谢 进程 密切 相关 。 
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3.2 HT BETEZE FEICOEDRLUR 

基于 KEGG Pathway 数据 库 对 差异 代谢 物 富 集 分 析 ， 注 释 到 68 条 代谢 通路 ， 其 中 3 
条 极 显著 的 代谢 通路 〈P< 0.01) 分 别 是 葵 丙 素 类 生物 合成 、 植 物 激素 的 生物 合成 和 类 黄酮 
生物 合成 。 其 中 茶 丙 素 类 生物 合成 是 最 显著 的 代谢 通路 , 茶 丙 氨 酸 为 下 游 次 生 代 谢 产物 提供 
原料 。 本 研究 基于 茶 丙 素 类 、 黄 酮 类 等 有 效 成 分 生物 合成 通路 共 往 选 获得 10 种 有 效 成 分 生 
物 合 成 有 关 的 差异 代谢 物 。 进一步 研究 发 现 , 茶 丙 素 类 生物 合成 相关 的 代谢 物 表 达 量 同 马 绢 
杜 鹏 花 的 其 他 时 期 相 比 , 在 衰老 期 和 凋谢 期 代谢 物 表 达 量 逐渐 升 高 , 而 黄酮 类 生物 合成 相关 
的 代谢 物 随 着 花 的 生长 发 育 含 量 旦 现下 调 表达 , 这 些 代谢 物 可 能 与 马 继 杜 鹏 花 凋 谢 有 重要 
控 作用 。 在 马 组 杜 鹏 花 的 花苞 期 茶 丙 氨 酸 表达 量 相对 较 低 ， 随 着 花 的 生长 发 育 代谢 物 表 达 量 
逐渐 升 高 , 代谢 通路 上 相关 代谢 也 发 生 相 应 的 变化 , 说明 它 对 整个 代谢 通路 的 完成 起 到 重要 
的 调控 作用 。 
现 有 研究 发 现 ， 葵 丙 素 类 生物 合成 是 最 重要 的 植物 次 生 代 谢 途 径 之 一 (Dong & Lin, 
2021) ， 植 物 的 茶 丙 素 类 生物 合成 谢 途 径 及 其 分 文 途径 会 产生 许多 次 生 代 谢 产 物 ， 如 色素 、 
酚 酸 、 黄 酮 、 类 黄酮 、 木 质 素 等 ( 文 欢 等 ，2017) 。 款 冬 花 不 同 生长 发 育 阶段 的 次 生 代谢 
| 物 代谢 组 成 明显 不 同 , 茶 丙 素 类 代谢 物 在 发 育 初期 至 中 后 期 表达 量 较 高 ; 黄酮 类 代谢 物 在 发 
~ 育 的 各 个 阶段 表达 量 均 有 波动 , 但 总 体 变 化 不 大 ， 次 生 代 谢 物 累积 促进 花 的 发 育 ( 贾 岩 等 ， 
i 2017) 。 采 用 高 通 量 测序 技术 对 川 佛手 不 同 发 育 时 期 转录 组 学 研究 ， 发 现 葵 丙 素 生物 合成 中 
7 个 相关 基因 随 着 果实 的 发 育 及 变色 呈现 显著 的 差异 表达 , 提高 其 表达 从 而 促进 果实 的 发 育 
QER 等 ，2020) 。 茶 两 素 类 生物 合成 途径 的 前 端 茶 丙 氨 酸 转化 为 反 式 肉桂 酸 ， 再 由 肉桂 
酸 羟 化 酶 〈C4H) 产生 对 状 基 香 豆 酸 ， 再 经 4- 香 豆 酰 -CoA 连接 酶 (4CL) 催化 生成 4- 香 豆 
酸 CoA， 进 一 步 向 下 游 特异 性 合成 途径 ， 转 化 为 不 同 的 葵 丙 素 类 代谢 产物 ， 包 括 香 豆 素 、 
黄酮 类 、 艾 类 和 木质 素 等 ( 王 玉 等 ，2019) 。 在 马 继 杜 鹏 花 开 花 至 凋谢 进程 中 葵 丙 素 类 生 
物 合成 相关 的 差异 代谢 物 苯 丙 氨 酸 (L-phenylalanine) 、 反 式 肉 桂 酸 Ctrans-cinnamic acid) ~ 
对 香 豆 酰 莽 草 酸 (p-coumaroyl shikimic acid) 、 松 柏 醇 〈coniferyl alcohol) 、 芥 子 酸 Csinapic 
acid) 、 紫 丁香 苷 〈syringin) 代谢 物 随 着 花 的 生长 发 育 含 量 呈 现 上 调 表达 ， 特 别 是 在 凋谢 期 
(F) 达到 最 高 值 ; 而 黄酮 类 生物 合成 相关 的 代谢 物 柚 皮 素 Cnaringenin)、 榭 皮 素 Cquercetin) 
和 查 尔 酮 (chalcone) 在 花 的 生长 发 育 过 程 中 下 调 表达 。 对 不 同 发 育 时 期 大 白 荣 花 的 转录 组 
研究 ， 发 现 花香 的 主要 成 分 是 葵 丙 素 ， 而 黄酮 类 化 合 物 是 花 、 果 实 和 种 子 的 主要 显 色 物 质 ， 
此 外 它 还 能 影响 花瓣 的 伸 长 ， 促 进 花粉 的 萌发 〈 李 改 珍 等 ，2017) 。 在 茶 丙 素 类 生物 合成 
ELANA REF 等 ，2015) 、 柚 皮 素 CTuetal,20160 、 阿 魏 酸 (Louie etal., 2010) 
等 作为 中 间 体 , 它们 的 表达 量 高 低 直 接 限制 到 下 游 代谢 途径 的 进行 ,从 而 影响 花 的 开花 至 调 
谢 过 程 。 葵 丙 素 类 、 黄 酮 类 等 代谢 物 作 为 重要 次 生 代 谢 产 物 ， 推 测 其 可 能 为 马 继 杜 人 鹏 花 开 花 
至 凋谢 进程 有 重要 的 调控 作用 。 为 了 更 好 地 了 解 代 谢 物 的 生物 合成 与 马 继 杜 鹏 花 开 花 至 凋谢 

进程 中 的 关系 ， 还 需 通 过 进一步 的 深入 分 子 水 平 研究 。 
本 研究 通过 非 识 向 代谢 组 学 技术 对 马 绢 杜 鹏 花 开花 至 凋谢 进程 中 不 同 阶 段 的 代谢 成 分 
进行 分 析 ， 基 于 KEGG 数据 库 对 差异 代谢 物 进 行 注释 并 获得 了 相应 通路 ， 构 建 了 马 缕 杜 鹏 
花 中 茶 丙 素 类 、 黄 酮 类 等 有 效 成 分 生物 合成 途径 ,为 进一步 研究 代谢 物 的 基因 调控 作用 以 及 
全 面 解析 代谢 途径 提供 代 参 考 , 为 延长 马 继 杜 鹏 花花 期 调控 机 制 提供 了 理论 基础 和 科学 依据 。 
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